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Abstract 

The diffractive process 7P — > V + X (where F is a vector me- 
son and X results from the dissociation of the proton) is studied. In 
particular, we consider the region of large momentum transfer (i.e. 
\t\ » ^qcd) large centre-of-mass (CM) energy, s. The asymp- 
totic (s oo, s/\t\ ^ 1) behaviour is derived from the BFKL equation 
and compared to that which is obtained in the Born approximation 
(two-gluon exchange). We also calculate the corrections to the Born 
graphs by iterating the BFKL kernel numerically. Improved conver- 
gence of the BFKL series is found by summing the logarithms which 
occur when an exchanged gluon goes nearly on shell. Importantly, 
we hnd evidence that the asymptotic solution to the BFKL equation 
is inappropriate over most of the HERA range and we provide more 
realistic predictions for the cross section. The predicted cross section 
is not too small and can be measured at HERA, up to momentum 
transfers Itl ~ lOGeV^. 



^Work supported in part by the Russian Fund of Fundamental Research. (93-02-03145) 
and by the Volkswagen- Stiftung. 



1 Introduction 



We consider the diffractive production of vector mesons in 7p interactions, 
where the momentum transfer —t ^ ^qcd the photon can be either 
real or virtuaL The vector meson (V^) is produced quasi-elastically while the 
proton dissociates into a diffracted mass, M (this contribution dominates over 
the elastic channel at large t). The scattering is mediated by the exchange 
of a colour singlet object, i.e. the perturbative ('hard') pomeron [|l], The 
largeness of the momentum transfer is expected to guarantee the dominance 
of short distance physics and hence the use of perturbative QCD. As we shall 
show, the reaction can be measured up to —t ~ 10 GeV^ at HERA and so 
we have the possibility to study the hard pomeron in detail. In particular, 
there is the prospect to examine the QCD pomeron trajectory away from 
t = 0. Elastic (i.e. without proton dissociation) vector meson production by 
photons at low momentum transfers has been studied in refs.|P, ||, ||. 

We start by considering the amplitude corresponding to Feynman graphs 
like those of figs.l(a,b). The contributions from the more complicated graphs, 
such as those in figs.l(c,d), will be discussed in section 4. Thus we can 
factorise the cross section into a product of the usual parton distribution 
functions and a quasi-elastic 'hard scattering' cross section: 

(1) 



and 



dt V TT / (iV2-l) 



d^kd^k'f'^{k,k',y)\ (2) 



Here, and throughout the paper, momentum vectors are transverse mo- 
mentum two-vectors of positive norm. We follow Mueller and Tang in 
choosing the normalisation such that if f^{k, k', y) is the amplitude for scat- 
tering two gluons (transverse momenta k and k') then the cross section of 
eq.(2) is that for quark-quark elastic scattering. The momentum transfer 
—t = Q^. The separation in rapidity between the final state parton and 
meson is y = ln(s/4g^) with s = x's {s is the 7*p CM energy). We define 
q^ = q\\+ Q"^ and = (Q^ + rriy)/^ {—Q^ is the photon virtuality). We 
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use —t and interchangeably throughout the paper but it should not cause 
confusion. 

For vector meson production, the amplitude is the product of the two- 
gluon scattering amplitude (see eq.(5)) and the form factor associated with 
the 7*V^ vertex: 



1 



2g2 2g2 + 2{k - Q/2)2 



(3) 



where 



31! 



Op 



This is the result of ref. (which is valid for heavy vector mesons, e.g. 
J/\I' and T). It is written in terms of the electronic width (r^g_) of the 
meson, so as to cancel a large part of the 0(q;s )-corrections. As usual, aem = 
1/137, Cp = 4/3 and Ca = Nc = 3. Although eq.(l) is written for heavy 
photon-proton interactions, it is also valid in the photoproduction limit, i.e. 

= 0. Also, eq.(l) is the component of the cross section corresponding 
to the scattering of transversely polarised photons only. The longitudinal 
contribution is obtained using cr^ = {Q"^/ My)aT and so it dominates at large 
enough Q"^. Recall that to get the electroproduction (ep) cross section one 
has to multiply eq.(l) by the photon flux, i.e. 



dx dQ 



IT X 



l-Y + Y'^/2) 



where x is Bjorken-x and Y is the photon energy fraction (= Q'^/{AxEeEp)). 

For the case of quark-quark elastic scattering, the form factor of eq.(3) is 
replaced by unity and the asymptotic BFKL behaviour of the cross section 
has been calculated to be 



da{qq qq) 
dt 



exp {2ujQy) 



4t2 [la^CAmvY 



(4) 



Where C(3) ~ 1.202 is the Riemann zeta function and ujq = In 2. In 

the next section, we calculate the cross section for diffractive vector meson 
production in the same (asymptotic) limit. We obtain analytic expressions in 
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the hmits \t\ ^ and \t\ qf- We compare the results with those obtained 
in the Born approximation (i.e. the two-gluon exchange model of Low and 
Nussinov [Q). These asymptotic expressions are only valid for ^ 1, i.e. 
the rapidity interval should be larger than the inverse QCD coupling. For 
smaller y, we expect significant deviations from the asymptotic formulae, 
and in section 3 we look for such deviations. Specifically, we iterate the 
BFKL kernel numerically, so as to obtain corrections to the Born term and 
to investigate the convergence of the BFKL series (recall it is an expansion 
in ~ asv). We perform our calculations up to and including terms ~ (osl/)^- 
We also show how to improve the convergence of the series by performing 
an all orders summation of the logarithms which are large in the asymmetric 
configurations where one (or more) of the t-channel gluons is nearly on shell. 
This resummation allows us to make predictions for the diffractive cross 
section for larger values of y than if we had used the simplest expansion. Of 
course, for sufficiently large y, even the resummed series is poorly convergent 
and a full solution to the BFKL equation is required. 

It should be acknowledged that, contrary to the case of elastic hadron- 
hadron scattering 0, the lowest-order photon-hadron amplitude is not sen- 
sitive to Sudakov corrections. It was shown in ref.[l^ that the Sudakov 
logarithms, which arise when internal quarks or gluons go on-shell, do not 
occur due to the pointlike nature of the photon. In addition, due to the 
largeness of g^, we need not worry about vector meson dominance contribu- 
tions to the vector meson production vertex. This is the case even in the 
photoproduction limit of = 0. 

In section 4 we return to the diagrams of figs.l(c,d). Such diagrams 
should be important at large M^(~ Q'^{l/x' — 1)). A summary of our results 
and predictions for cross sections in the HERA range are presented in the 
conclusion. 



2 High energy asymptotics 

The solution for the QCD pomeron amplitude at non-zero momentum trans- 
fer was given by Lipatov in terms of the eigenfunctions in the transverse 
coordinate (p) representation, i.e. 

f'ik, k', y) = T^J dy .,l\i^.2 Q)I^{k, Q), (5) 
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where 



P1P2 



[ik-pi+i(Q-k)-p2] 



(6) 



The co-ordinate vectors of the t-channel gluons are pi and p2 (they are the 
vectors conjugate to the transverse momentum vectors k and k — Q) and 
V{k, Q) is the 'impact factor' for the vertex, A. For the couphng to a quark 
hne, one can take an impact factor of unity with the modification of eq.(6) 
which is appropriate for scattering off coloured particles (this is the prescrip- 
tion of ref. [0] and we shall discuss it in more detail in section 4). Whilst for 
the 7*1^ vertex we take the form factor of eq.(3). The function uj{i') is 



2a,CA 



TT 



7^e[7/;(l) - V^(l/2 + iz/)] 



with ipi^) = {d/dz)\nr{z) (the derivative of the logarithm of the gamma 
function). 

In the high energy limit y —>■ 00 and the dominant contribution comes 
from the saddle point at i/ = 0. Thus we are interested in the function Iq. 
In the qq case it was calculated in ref. 



-^I'oik,Q) - ^. 



(7) 



In the 'j*V case, things are a little more complicated. Let us define 



V{p) 



d^k 
(2^ 



J{k-Q/2)-p 



2g2 2g2 + 2{k - Q/2y 



(8) 



R = {pi + P2)/2, p = Pi — p2 and Kq is the McDonald function. Now instead 
of eq.(7), we have 



/^^(g) = C j d^Rd^p- 



R + p/2\\R- p/2\ 



JQ-R 



V{p) 



-C J JoiQR) RdR 



\P\ 



\R + p/2\\R~ p/2\ 2 



(9) 
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and Jo{x) is the Bessel function of the first kind. The cross section is then 
given by 



X 



47r \ 



Q 



(10) 



where the quark-quark scattering cross section is given in eq.(4). Eq.(9) can 
be evaluated analyticallyf] in the limits ^ and <^ gy. 
ForQ^^gJ (see ref.|lTl) 

C_ ( \ 327r- 



/o^(Q)b^»,^--7%ln|^lT;^^. (11) 



In the other limit of small <^ gy , the essential contribution comes 
from R ^ p/2. The logarithmic integral over R goes from p/2 to 1/Q and 
is equal to ^^12 ^ ~ In while the integral d'^pp Ko{pq\\) = vr^/g|. So 

/0^(Q)lQ^«.t--|3l^(^)^'- (12) 

In fig. (2), the results of an exact numerical computation of Iq (eq.(9)) are 
compared with the asymptotic formulae of eqs.(ll) and (12), as a function 
of the ratio r = (5^/4g^. This is a sensible variable to plot all our results 
against, since it exploits the corresponding scale invariance of the BFKL 
kernel. 

To finish the section let us compare the asymptotic BFKL behaviour for 
the QCD pomeron amplitude with that obtained in the Born approximation. 
For the two-gluon pomeron the cross section for 7*g — > Vq takes the form 

^"'^^^^ ^ = c^i^ J^ (13) 

dt 81 ' ^ ^ 



^We are indebted to Hans Lotter for correcting a mistake in our initial evaluation of 
eq.(ll) 
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where 



J 



TT J 



(Pk 



k'^iQ - kf 



1 



2^2 2(i + 2{k - Q/2f 



- 8 1 + .... 

For large ^ the leading logarithmic contribution to J comes from 
-C (A; - < (5^/4 and J7 = ^In The same logarithmic be- 

haviour arises in the asymptotic BFKL calculation (eq.(ll)), which contains 
an additional enhancement due to the factor ~ exp (2a;o|/)/|/^. Note that 
at not too large energies, the numerical coefficient of the BFKL asymptotic 
result can be small and for uj^y 1.5 it is only about 3% of the Born term. 
For small < gy the BFKL factor of 



Q 

is to be compared with the 



In^ exp {2u^y)l[-asCAam 



1 , 2 



obtained in the Born approximation. There is an additional enhancement 
of ~ Q^/Q'^ in the case of the BFKL pomeron. The origin of the difference 
is clear. The distance between two gluons in the Low-Nussinov pomeron is 
fixed since they must couple to the qq pair, i.e. their size docs not exceed that 
of the upper vertex (p ~ V^ll)- Thus the two-gluon system is not sensitive to 
small Q^. However, in the BFKL pomeron, after a few iterations (or rungs 
of the ladder) the gluons can be separated by large distances p ~ l/Q and 
so the amphtude ~ l/Q. 



3 Iterations of BFKL kernel 

So far, we have considered only the asymptotic solution to the BFKL equa- 
tion. In this section we consider ordcr-by-order iteration of the BFKL kernel 
for non-zero t and examine the nature of the leading log s expansion. For 
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sufficiently small z (= ijCac^s/'^'^) (and sufficiently large \t\) we can expect 
the Born approximation to be appropriate. As z increases, so the need to 
include more and more terms in the BFKL series increases until the point 
is reached where an all orders summation is vital. By studying the BFKL 
series in this way we can hope to improve on the Born cross section estimate 
in the region of intermediate z. In addition, such an expansion will allow us 
to identify the onset of the region where the BFKL (all orders) resummation 
is vital. To start, let us recall the BFKL equation for non-zero t: 




For iterations from the quark line we use the boundary condition: 

Mk,z) = l. (16) 

As usual, all the hi are two- vectors in the transverse plane and k — k — Q. 
The first iteration of this input leads to the simple expression: 

Mk,z) = zlnl^-^y (17) 

Subsequent iterations are difficult to perform analytically and we evaluate 
them numerically. 

Successive iterations lead to large oscillatory behaviour in the regions 
/c^ <^ and /c^ ^ (as a result of logarithms like the one in eq.(17)). 
This behaviour leads to a very poorly convergent series. To overcome this dif- 
ficulty we will sum up the (double) logarithmic contributions ~ [2;ln(/c^/Q^)]" 
analytically, and then factorise these badly oscillating terms to leave behind 
a better convergent series. 

There arc no infrared divergences in eq.(15) (neither at k' — > k, nor at 
k' or k' Q) and the only logarithm comes from the last (reggeization) 
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term in the region k"^ <^ k'"^ <^ . For small k'^ ^ (or k"^ <^ Q"^) eq.(15) 
leads to 



d^{k,z) 
dz 



£k' 



nJ {k'-kf 
-$(A;,z) 



{k' - kf 



k'^ 



+ 1 + 



..l2 



$(fc',^) 



+ 1 



_k'^ + {k' -kf 

and we have introduced the function $(A;, 2;), where 

n=oo 

^{k,z)= J2 ^n{k,z). 



(18) 



n=0 



We have used the result that 



0. After the angular integration 



d^{k,z) 
dz 



dk'- 



^k,z) 



+ 

1 + k^k'^ 



^k',z) 



kyk'^ 



Which can be written more concisely: 



{^{k',z) -^{k,z)) + ^{k,z) 



(19) 



(9$(fc, z) 


ff{ 




dz 




[le-ei J 



1 



(20) 

where ^ = k^/Q'^ and ^' = fc'^/Q^. It is easy to see that the only logarith- 
mically large contribution comes from the last term, which givesQ ^{k, z) In^ 
and so 



d^{k,z) 
dz 



i d^'mk',z) -^k, z)) 



•^At first sight the accuracy of eq.(18) is insufficient to be sure we are not missing some 
constant in comparison with the In^. However this log comes only from the reggeization 
part of the BFKL kernel, for which the exact answer is known. The contribution is equal 
to z)[ln((5^//i^) + ln(fc^//x^)]. The infrared cutoff is cancelled after we subtract 
the part 



2$(fc,z) 







Jo 





= 2$(fc,z)ln(fcV^2)^ 



which has already been included in the integral of eq.(21). Finally one gets $(fc, z)ln^ 
and the accuracy of eq.(20) is of the order of 0{k^ /Q^). 
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In the double log approximation the solution is 



^{k, z) = exp {zln^) 



(the integral over ^' does not give rise to any logs). Now we can put $(A;, z) = 
(p{z,^) exp (zln^) and compute the solution for (p 



with the initial condition 0(0, = 1 (corresponding to iterations from the 
quark line). However one can only make sense of such a framework for z ^ 1. 
At larger z the last integral in eq.(22) diverges in the region of ^' — > oo. The 
implication is that, for z ^ 1, eqs. (18-22) are not self-consistent, i.e. the main 
contribution to $(A;, z) (with small ^) comes from the region of large The 
results of the numerical calculation are presented in fig.3. For z ^ 0.4 the 
variations of the function (p (from unity) are rather small. 

Let us explain how the diffractive cross section is obtained from the 
^m{k,z) functions. First define (and we now make explicit the dependence 
upon the scaling variable r = Q'^/Aq?,) 



where il)^{k)z'^/m\ — $^(A;,2;) is a dimensionless function corresponding 
to m iterations of the BFKL kernel starting from the quark line boundary 
condition, ipQ^k) = 1. Iterations from the vector meson boundary condition 
of eq.(3) lead to the functions ip^{k). The Born term is just the lowest order 
convolution of eq.(14). The Born term and first iteration from the quark 
line are known analytically, and the second iteration has been computed 
numerically. Similarly, the zeroth iteration from the vector meson boundary 
condition is known analytically, whilst the first iteration can be computed 
numerically etc.. Consequently, we can convolute these results to obtain 



dz 



2 





(22) 




(23) 
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corrections to the Born contribution. Where possible we have checked our 
results by calculating the convolutions of eq.(23) in different ways, i.e. we 
have checked that ipf ® ipY = '4^2 ^ "^o- For the cross section we have 



da'^{Yq Vq) Ana 



(24) 



dt 81 
After summing the double logarithms 

oo n oo n 

and the equality means that 

<i>m = rn{k)-rn-i^nii--—+ +v;L^(fc)in-a / \: , + 

[n — 1)\ [n — mjlml 

We choose to resum terms ~ -zln^^ (this expression simultaneously reduces 
to the double log solution when ^ <^ 1 or ^ ^ 1). By convoluting the 
(pl functions (rather than the ip!^) with ipQ, in eq.(23), we can compute the 
double log improved cross section. It should be appreciated that it is only 
meaningful to perform the exponentiation in the small region. Since the 
integrals (in the convolutions) are performed over all ^ (up to infinity), we 
must introduce a 'factorisation' scale, to delineate the small region, 
i.e. < ^0 must be satisfied before resumming. Choosing = therefore 
corresponds to the un-resummed series. We choose = 1 but the choice 
is essentially arbitrary (in the same way that the QCD factorisation scale 
is arbitrary). Essentially, the larger one chooses C,o, the better the effect of 
resummation - up to the point where the logarithm of can no longer be 
considered 'large'. 

In figures 4 to 6 we show the dependence of the function ip{z, r) upon r for 
different values of z (in fact we plot [ip{z,T)t'^ / {20)]^). These are essentially 
plots of the cross section, which can be obtained from the plots via eq.(24) 
and eq.(l). Notice that we have removed the explicit 1/t^ behaviour of the 
cross section to leave behind a quantity which only depends upon the 'scaling' 
variable, r, and the 'energy' variable, z. In all cases, the dotted lines are the 
BornQ results, the dashed lines include the 1st (~ z) BFKL corrections (to 

'*0r the Born plus the double logs in the case of fig. (5,6). 
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the Born results), the sohd hnes include the 2nd (~ z"^) BFKL corrections 
and the diamonds include the 3rd (~ z^) BFKL corrections. The dash-dotted 
lines are obtained using the asymptotic results of the previous section (i.e. the 
numerical evaluation of eq.(9)). In fig. (4), we show results for the simplest 
expansion of the BFKL series, i.e. without the double log resummation. 
Fig. 4 (a) shows quite clearly the good convergence that is expected for small 
enough z {z = 0.05). The Born term provides a reasonable description 
and lies within 30% of the higher order result over most of the r range. 
The major deviations arise around the dip near r = 1 (the cross section is 
zero here in the Born approximation due to the cancellation between the 
diagrams corresponding to the two gluons coupling to the same and different 
quark lines). Going to z = 0.35, as we do in fig.4(b), we find very poor 
convergence for r > 1 (where the 2nd and 3rd iterations often deviate by an 
order of magnitude). The problem is a result of the oscillatory nature of the 
solutions to the BFKL equation. The first corrections to the Born term are 
big and they completely fill in the dip. The next corrections tend to restore 
the dip, by including large negative contributions in the region of r > 1 which 
partially cancel the positive contribution of the first corrections. This leads 
to the unpredictable behaviour of fig. 4(b). We do not expect, therefore, to 
make reliable predictions for the diffractive cross section using such a simple 
minded expansion of the BFKL series (for all but the smallest values of z). 
Before moving on to the results obtained with the resummed series we should 
point out that, as noted above, the poorest convergence of the un-resummed 
series occurs for the larger values of r, (in particular in the region around 
the dip) and that this is a general feature of our results. Indeed we shall 
see that it is precisely this region that is improved by resumming. This is to 
be understood as the region where the asymmetric configuration (i.e. where 
one of the gluons which couples to the ^yV vertex carries all the momentum 
transfer) is important. 

The double log improved results are shown next. Fig. 5 shows the anal- 
ogous plots to those in fig.4. The accuracy of the (resummed) Born term is 
improved ai z = 0.05, whilst at z = 0.35 the higher corrections are necessary. 
However, for z — 0.35, we can now claim to make reasonable predictions over 
the whole r range after just two iterations of the BFKL kernel, i.e. the dif- 
ference between the 2nd and 3rd iterations is typically < 20%. Fig. 6 shows 
the resummed results for z = 0.2,0.5 and 0.8 (which are typical values for 
HERA). For z = 0.2 we obtain convergence at the level of ~ 1% at order 
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. Going to larger 2;, as one expects, the quality of convergence deteriorates 
and, at z = 0.8, we clearly need a more complete summation of the BFKL 
series. However, it is essential to realise (since the cross section is expected 
to vary rapidly with the momentum transfer) that a prediction of the cross 
section to within a factor ~ 2 is very useful. For example, at z = 0.8, the 
Born cross section lies typically more than an order of magnitude (~ 50) be- 
low the order z^ results (diamonds), which is in turn a factor ~ 2 — 3 below 
the asymptotic results. 

Comparison of the numerical results obtained in this section with the 
asymptotic results (dash-dot lines) shows quite clearly that the asymptotic 
formulae do not represent the diffractive cross section over most, if not all, 
of the kinematic range examined. Only at the largest value of z, does it 
look likely that, after fully summing the BFKL series, the result will agree 
with the asymptotic prediction. In the region where the double logarith- 
mic resummation is most important, i.e. r > 1, the asymptotic formulae fail 
completely to approximate even the gross features of the expected shape, 
e.g. the dip. In addition, the asymptotic prediction always appears to over- 
estimate the size of the cross section. Ultimately, for large enough z, the 
asymptotic result must be more appropriate than the 'fixed order' results. 
We can see, in fig. 6, that at z = 0.8 the BFKL series is slowly convergent and 
one can no longer trust the order z^ predictions, i.e. this marks the onset 
of the dynamics manifest in the full BFKL solution which will ultimately be 
described by the asymptotic solution. It should not come as a surprise that 
the asymptotic prediction is not trustworthy for z ^ 1. Its reliability relies 
upon the validity of the approximation that the exponential in eq. (5) can be 
expanded about u = 0. This is a very poor approximation for z ^ 1 since 
the integrand of eq.(5) actually vanishes at z/ = and rises rapidly to its 
maximum value only for large z. 

It is also worth commenting upon the results of Lipatov and of Hancock 
and Ross |T^. They have found that, in the case of a running as{Max{k'^, Q^)), 
the leading eigenvalue of the BFKL kernel lies somewhat below its asymptotic 
value for all realistic energies (and non-zero t). Therefore the cross section 
is expected be smaller than that predicted by the asymptotic solution based 
on eq.(4). We do not investigate the effect of a running coupling here. 
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4 Sub-leading Corrections 



We have devoted a whole section to the calculation of the process 7*g — > Vq. 

We now turn to a discussion of the 'lower part' of the diagram and the validity 
of the simple assumption that the pomeron couples to a single parton inside 
the proton. 

If is large enough, so that the partons (which we label by their lon- 
gitudinal momentum fractions x' and x") within the proton are separated 
by distances p' S> 1/Q then we are certainly entitled to restrict ourselves to 
the simple coupling to a single parton line. This may seem somewhat supris- 
ing, since at first sight eq.(6) appears to suggest that the only contribution 
arises from the coupling to different parton lines (the amplitude vanishes at 
Pi — P2)' This is not the case. To exploit the conformal invariance, Lipa- 
tov redefined the eigenfuctions (by adding terms proportional to S^'^^k) and 
5'^'^\k — q) which vanish on coupling to colourless particles) and so the con- 
tribution from the Feynman graphs where the pomeron couples to a single 
parton is hidden. Indeed, let us start by assuming that p' ^ such that 
(for large Q^) 

J \R + p/2\\R- p/2\ )^ \ ) 

i.e. the main contributions arise from the regions where R is within A/? ~ 
l/Q of the singular points at ±p/2. This is the result which underpins 
the Mueller- Tang prescription for calculating the contribution from the cou- 
phng to a single parton line. The first term of the r.h.s. of eq.(25) corre- 
sponds to interaction with a quark at p' = p/2 whilst the second term is 
due to interaction with an anti-quark at —p/2. Corrections to eq.(25) are 
~ (Ai?/p)^ ~ l/(p^(5^) and represent contributions from Feynman graphs 
where the pomeron couples to more than one parton. 

Note that this approach is not appropriate for the ^yV vertex. In evaluat- 
ing the R integral of eq.(9) we are not entitled to assume that the essential 
contributions arise from the regions around the singular points since the sub- 
sequent p integral would be dominated by contributions from small p ~ 
i.e. after evaluating the it! integral using eq.(25), we would be left with 

/ Ko(q\\p)e"^-'/'£p 
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where the dominant contribution arises from small p ~ 1/Q and thus we 
are not justified in assuming that p ^ Ai?. The corrections are therefore 
essential and indeed give the main contribution to Iq{Q) (i.e. the logarithm 
in eq.(ll)) at large Q"^. 

At the lower vertex, the situation is quite different since the target disso- 
ciates. The impact factor now takes the form, 



|/P(Q)|2 = 2 j d^p d^R d^E! |^(p)p 



p e^<5-« p e-^«-^' 



\R- p/2\\R+ p/2\ \R - p/2\\R' + p/2\ 

(26) 

where |\l/(p)p is analogous to V{p) in eq.(8). The integrals over R and i?' can 
now be performed safely in the p ^ Ai? limit since the exponentials {e^^'P/"^ 
and q~^Q'pI'^') cancel each other. Thus we get 

\im?=^-^jd'p\^{p)?- (27) 

With the lower limit of the p integral equal to 1/Q and the upper limit fixed 
by the proton radius the integral over \1/ gives the conventional structure 
function, i.e. ~ 1/p^ generates the leading logarithmic behaviour of 

the parton density function with the typical values of p ^ ^IQ- Since the 
integrand of the p integral is independent of (it is simply the wavefunction) 
it follows that the main term in eq.(26) is indeed the contribution from the 
graphs where the pomeron couples to a single parton line. Thus we have 
obtained the Mueller- Tang result that was used in the first part of the present 
paper and have justified our assumption that the dominant contribution to 
the proton dissociation (at large enough Q^) arises from the coupling to 
a single parton line. For a more detailed discussion of the nature of the 



pomeron-quark coupling we refer to ref. |T3 . 

Note that we neglected the cross terms ~ e*'^''^ which arise from the 
points R ±p/2; R' — » ±p/2 after the R and R' integrals of eq.(26). 
They represent contributions from the interference between graphs where 
the pomeron couples to single but different parton lines. Such contributions 
are suppressed at least by logarithms of since the momentum transfer 
must flow into the proton wavefunction. 

Similarly, diagrams like the one of flg.l(c) are beyond the leading InQ^ 
approximation since they lose at least one power of In due to the hooking 
of the large transverse momentum, /c ~ Q/2 into lower rungs of the 'structure 
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function'. Q Therefore the contributions from such graphs are suppressed in 
comparison with those of figs.l(a,b) by the factor — ^72- The anomalous 
dimension of the twist-2 operator (structure function) is 72 oc and its 
presence reflects the fact that one loses one power of InQ^. The additional 
smallness (—1/2) comes from the colour coefficient (C^/2 in comparison with 
Ca). Note that the coefficient is negative and so the graph in fig. 1(c) describes 
a colour screening effect. If the BFKL gluon continues down and touches the 
next s-channel gluon (i.e. x'") then one gets a factor of 1/4 and two powers 
of 72, and so on. When computing the cross section, one must remember to 
double the contributions from configurations where the two t-channel gluons 
couple to different s-channel gluons. Finally, if the two gluons couple to the 
same s-channel gluon {x") then one loses one logarithm and picks up a colour 
coefficient which is the square of the result obtained when the gluons couple 
to different gluons. Thus the leading contribution to the order corrections 
is suppressed relative to the contributions of figs.l(a,b) by the factor 

7^2 = (-^ - i + I) 72 = -^72 ^ -0.3. (28) 

We put as = 0.16 and use the double logarithmic approximation result: 
72 = Cao^s/ (t^^) — 0.36 with u = ujq = 0.42. Such contributions correspond 
to the ~ l/(p^(5^) corrections to the impact factor, Iq{Q), calculated above. 

Let us now turn to graphs like the one shown in fig. 1(d). Such diagrams 
differ from those of fig. 1(c) in the sense that one of the gluons now couples 
deep inside the proton (e.g. to a different branch of the parton cascade). As 
such we expect a power like suppression combined with an enhancement due 
to the sampling of the two-parton component of the hadron wavefunction. 
Let us estimate the possible size of such contributions. The contribution is 
analogous to the graph of fig. 1(b), in the sense that the pomeron is coupling 
to different parton lines in the scattered object. Unlike the case of fig. 1(b) 
however, its contribution is not calculable in perturbation theory since it is 
dependent upon the rather complicated structure of the proton wavefunction. 
Since the proton is more diffuse than the compact quark-antiquark pair (of 

^It was shown in section 2 that for large y the internal momenta inside the BFKL 
ladder are typically k Q/2. The contribution from the region of small k (or [Q — k)) 
dies out with energy. Even in the Born (two-gluon) approximation the main logarithmic 
contribution comes from {k — Q/2) ^ Q/2 , i.e. k close to Q/2. 
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the upper vertex) we anticipate that fig. 1(d) will be suppressed at all values 
of the momentum transfer considered here, i.e. Q ^ Aqcd- 

In comparison with eq.(l), the contribution of the graph of fig. 1(d) to the 
cross section contains an additional colour factor 1/{N^ — 1) = 1/8. There is 
additional suppression from the integration over the transverse momenta of 
the t-channcl gluons, ki and k2{= Q — ki), which connect the pomcron and 
the proton. In the case of fig.l(a,b) there were no correlations between the 
momenta {ki , /C2) in the amplitude A and the corresponding momenta {k[ , 
k2) in the complex conjugated amplitude A*. One therefore has independent 
integrals over ki and k[ which give two large (of the order of Q^) factors. 
This is not the case for the diagram of fig. 1(d). The transverse momentum, 
k2, is typically carried away by its associated final state gluon which means 
that k2 ~ k2 and we get only one large integral which is ~ Q^. The remain- 
ing integral is then over the intrinsic gluon momentum, dk (ki^ — k2 + Sk) 
and we become sensitive to non-perturbative effects in the proton, i.e. the 
momentum transferred over the k2 loop in fig. 1(d) is limited by the wave 
function of the target nucleon. However, there is an enhancement due to 
the larger combinatorial factor associated with coupling the two-gluons to 
different parton lines. 

We can be more explicit by returning to the impact parameter represen- 
tation. As we have just shown, the square of the impact factor for coupling 
the pomeron to a single parton line is oc (47r/Q)^ (eq.(7,27)) and to obtain 
the cross section, we then multiplied by the parton density factor (which can 
be thought of as the multiplicity associated with the pomeron-proton impact 
factor). An analogous, but slightly more involved, calculation can be per- 
formed for the case of the diagrams like that in fig. 1(d), i.e. we can compute 
the square of the impact factor and multiply by the colour factor (1/8) and 
the square of the parton densities to obtain the cross section. The impact 
factor can be written thus: 



The Mueller- Tang subtraction term (denoted by 'M.T.') is defined by the 




M.T. . (29) 



R-p/2\\R + p/2\\R'-p/2\\R' + p/2\ 
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replacement: 

\R-p/2\\R + p/2\ J ^\R-p/2\\R + p/2\ Q^^ ' 

and similarly for the R' integral. In this way we remove the contribution cor- 
responding to the coupling to a single parton line. The answer, as it must be, 
is sensitive to the large distance physics of the nucleon and this is contained 
in the two-particle wavefunction '^2{p)- Dependence upon longitudinal mo- 
menta is assumed to factorise into the parton number densities. 

We use a gaussian form for the p distribution and can introduce the same 
scale, Qg, that is used by L.V. Gribov, Levin and Ryskin (GLR) to pa- 
rameterise the screening corrections to the deep inelastic structure functions 
SHI, i-e. 

|vl/2(p)|2 = ^e-'''«o/4 (30) 

with Qq = 1.2 GeV^ from fits to SppS data (mainly on the inclusive cross 
section for charged hadron production). 
Eq.(29) can be simplified, i.e. 

|AJo^(Q)P = 2j£p p'\<^2{p)? Up)Jy{p) (31) 

where 



lx(l — x) 







- r^, ^o(QPv^) (i + ,.Q.p). (32) 

Jl/2 \ X 



and similarly for Jy{p)- In the limit Q/Qo ^ 1 the integral is dominated by 
the contribution from x ~ l/{QpY <^ 1 and p ~ 1/Qq. We can thus neglect 
the integral over x > 1/2. After expanding the exponentials and performing 
the remaining integrals we obtain: 

\^m)? - '-^p (33) 
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with S ~ In^Q/Qo). We also computed S numerically, and for AQ^/Ql = 10 
we get S = 0.74. 

Combining eq.(33) with the colour factor and the square of the parton 
densities we can now estimate the size of the corrections which originate from 
fig. 1(d), i.e. we have, rather than the factor G{x' ,t)dx' , the term 



S fQl 



8 \4g2 



G{xi,Qy4:)G{x2,Qy4:)dxidx2. (34) 



Now let us talk a little more about the value of Qq. At first sight, one 
might anticipate that Qq ~ ^/Rn (the inverse proton radius). However 
from the semihard phenomenology of ref.[0 we know that the gluon-gluon 
correlation length Ro is much smaller than R^. Also, calculations based 
upon QCD sum rules determine the radius of the two-gluon form factor of the 
proton to be Rq ~ 0.3 — 0.35 fm |]T6[, which corresponds to Ql ~ 1.2 GeV^. 
We should point out that we have ignored corrections due to interactions 
between the two branches of the parton cascade (i.e. we used the square of 
the single parton density in eq.(34)). 

We are ready to make a numerical estimate for the size of this correction. 
The integrations over xi and X2, in eq.(34), are limited by the mass of the 
hadron system ~ -Ij- ( — + — ) or by other experimental conditions. For 
example, one might choose to search for jets in the proton dissociation. This 
would be a useful search to perform, since it would not only give an additional 
way to measure the momentum transfer (and the momentum fractions Xi , X2 
carried by the jets) but it could be used as a means to distinguish between the 
leading and sub-leading contributions (for the sub-leading contribution, the 
vector meson transverse momentum is balanced by a pair of jets each with 
transverse momentum ~ (5/2, while in the leading case it is balanced by the 
one jet). Unfortunately, it is not easy to make such a study at HERA since 
the need to ensure a large rapidity gap forces the jets to lie close to the beam 
hole. Since we do not require the observation of the proton dissociation, we 
shall estimate the sub-leading corrections assuming no jets are seen in the 
proton direction, i.e. 



e ~ ^^-^ < 3° 

xEp 
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in HERA lab frame then 



Xl,2 > Xm 




11.1 1 -J 



for \t\ = 3 GeVl 
tributions to be 



Thus we find the ratio of the leading to sub-leading con- 




(35) 



and we have put Xm = 0.02. Taking = 1.2 GeV^ and \t\ > 3 GeV^ (i.e. 
S fa 0.74) it follows that this ratio is < 5%. 

As we mentioned above, perhaps the best way to study the sub-leading 
contributions would be to look at the jet structure of the final state. At high 
enough centre of mass energies one can have a large diffracted mass and still 
maintain a rapidity gap. Events with a large diffracted mass will lead to 
relatively large values for the sub-leading corrections (since Xi and X2 can be 
smaU). 

Based on these estimates we think that, over the HERA range, the sub- 
leading corrections discussed here do not change the results of the previous 
calculations (which assume factorisation of the proton dissociation) in a cru- 
cial manner. The effect may be less than ^ 5% for dissociation into a small 
mass, but it would be interesting to study experimentally how it reveals itself 
at higher masses. 

5 Conclusion 

We have studied the hard diffractive production of heavy vector mesons in 
deep inelastic scattering. The hardness is provided by a large momentum 
transfer, \t\, and allows one to probe the essential dynamics which determine 
the pomeron of perturbative QCD. We have made predictions for the cross 
section for a range of kinematical configurations accessible at HERA. Cor- 
rections to the simple two-gluon picture of the pomeron have been computed 



19 



using the leading logarithmic formalism of BFKL and they are seen to be 
very large. 

It is encouraging that the total cross section for this process is not too 
small. Specifically, for ^/s ^ 200 GeV, photon virtuality ^ 10 GeV^, 
x' > 0.01 and |t| > 2 GeV^ (i.e. a mean z of 0.6) the total cross section 
a'^{'-f*p V -\-X) for the hard diffractive production of J/^ mesons is ~ 1.5 
nb. This is to be compared with the asymptotic BFKL prediction of ~ 3 nb 
and the Born prediction of 84 pb. Although this cross section (1.5 nb) is 
small, (it is only 0.01% of the total DIS cross section) it is not prohibitively 
so. 

Although we have concentrated on the general case of non-zero Q^, all 
of our results could equally well apply to photoproduction (where the pro- 
duction rate is much higher). We need not worry about Sudakov or VMD 
corrections due to the largeness of q'^. For example, for = 0, = 200 
GeV and \t\ > 2 GeV^ one is probing a mean z of 0.8 and the cross section 
for J/\l/ production off transverse photons is ~ 15 nb. At -^s = 100 GeV 
this cross section falls to ~ 5 nb. 

It is an important asset of this process that one does not need to observe 
the products of the proton dissociation in order to extract the momentum 
transfer, t. In DIS, observation of the scattered electron together with the 
decay products of the vector meson allows a clean determination of t whilst 
in photoproduction one can assume that the incoming photon is coUinear 
with the incoming electron and hence the px of the vector meson gives t. 
Consequently, the rapidity gap between the vector meson and the scattered 
parton can be as big as 9 units. As a result, HERA can probe the region 
where the BFKL series is poorly convergent, i.e. where the resummation of 
leading Ins terms is most important. 

We restricted ourselves to heavy mesons for reasons of simphcity although 
p (or (f)) production will also provide important information on the QCD 
pomeron with generally a much higher production rate. Our calculations 
have been totally general and so one could use them to predict the production 
cross sections for any vector meson. However, we expect the non-relativistic 
approximation to be less appropriate for the light mesons. In addition, VMD 
contributions may well start to become significant in the case of photopro- 
duction for not-too-large t. It is interesting to note that for \t\ ^ my the 
yield of heavy vector mesons is larger than that of the lighter ones. This 
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is because the cross section for scattering transversely polarised photons is 
proportional to T^+^-rriy (see eqs.(l),(3)). Since the electromagnetic widths 
of the J/\l/ and p are similar we therefore find the J/\i/ rate to be enhanced 
by a factor ~ 50 over the p production rate. This is a simple consequence 
of the fact that at large \t\, the pomeron has enhanced coupling to the more 
compact cc pair which constitutes the J/\l/. However, this enhancement is 
diluted to a single power of my in the case of longitudinal photon scattering, 
which dominates at large enough photon virtualities. We note that chang- 
ing the photon virtuality allows one to scan through the vector meson wave 
function (as was discussed in ref.[p!7|). 

Corrections to our predictions which arise from higher order QCD and 
non-perturbative effects (in the dissociation of the proton) have been es- 
timated and should not be essential over the HERA range, perhaps con- 
tributing by no more than ^ 5%. However, higher order corrections to the 
BFKL amplitude itself may be more significant. For example, if we take 
y = ln[s/(4Max(Q^, my, Q^))] (rather than y = ln(s/4g^)) then the cross 
sections quoted in these conclusions are reduced by ~ 40%. This ambiguity 
in the definition of y is a direct consequence of the leading log nature of the 
calculation. 

It should not be long before the HERA experiments are able to study 
data in the regime considered in this paper. In particular, there is already 
some preliminary data on both J/^f and p production out to |t| values as 
large as 1.5 GeV^ (in DIS and photoproduction) and at centre-of-mass 



energies of ^ 100 GeV |18 



Acknowledgements: We should like to thank Prof. M.Arneodo, Dr. H.Jung, 
Prof. C.Peroni and Dr. P.J.Sutton for many fruitful discussions. Very special 
thanks to Prof. J.Bartels, H. hotter and Dr. M.Wiisthoff for their interest 
in this work. One of us (M.R.) gratefully acknowledges the hospitality of 
the Torino University Physics Department, the Theory Group of the RAL 
Particle Physics Department and the DESY Theory Division, where part of 
this work was done. 



21 



References 

[1] V.S.Fadin, E.A.Kuraev, L.N.Lipatov: Sov.Phys.JETP 45 (1977) 199; 
la.Ia.Balitski, L.N.Lipatov: Sov.J.NucLPhys. 28 (1978) 822 

[2] L.N.Lipatov: Sov.Phys.JETP 63 (1986) 904 

[3] A.Donnachie, P.V.Landshoff: Phys.Lett. B185 (1987) 403 

[4] J.R.Cudell: NucLPhys. B336 (1990) 1 

[5] M.G.Ryskin: Z.Phys. C57 (1993) 89 

[6] S.J.Brodsky et al.: Phys.Rev. D50 (1994) 3134 

[7] A.H.Mueller, W.K.Tang: Phys.Lett. B284 (1992) 123 

[8] F.E.Low: Phys.Rev. D12 (1975) 163; S.Nussinov: Phys.Rev. Lett. 34 
(1975) 1286 

[9] J.Botts, G.Sterman: Phys.Lett. B224 (1989) 201; Nucl.Phys. B325 
(1989) 62 

[10] G.R.Farrar, G.Sterman, H.Zhang: Phys.Rev. Lett. 62 (1989) 2229 
[11] J.Bartels et al.: in preparation. 

[12] R.E.Hancock, D.A.Ross: Nucl.Phys. B394 (1993) 200 

[13] J.Bartels et al.: 'How does the BFKL pomeron couple to quarks?', 
DESY-94-244. 

[14] L.V.Gribov, E.M.Levin, M.G.Ryskin: Phys.Rep. 100 (1983) 1 

[15] E.M.Levin, M.G.Ryskin: Phys.Rep. 189 (1990) 267 

[16] V.M.Braun, R.Gornicki, L.Mankiewicz, A.Schafer: Phys.Lett. B302 
(1993) 291 

[17] B.Z.Kopeliovich, J.Nemchick, N.N.Nikolaev, B.G.Zakharov: Phys.Lett. 
B309 (1993) 179 



22 



[18] ZEUS CoUab., M. Derrick et al.: 'Observation of vector meson produc- 
tion in electron-proton DIS at HERA', ICHEP94 Ref. 0663; ibid. 'Cross 
Section Measurement of the process ep epJ/'^ at HERA', ICHEP94 
Ref. 0672; I.Abt et al.: HI Collaboration, 'Photoproduction of J /^^ 
Mesons at HERA', DESY 94-153. 



23 



Figure Captions 



Fig. 1 Feynman graphs corresponding to diffractive vector meson produc- 
tion in DIS. 

Fig. 2 The asymptotic prediction for the diffractive scattering amphtude. 
The dotted hnes are the analytical approximations and the solid line is 
the result of exact numerical computation. 

Fig. 3 The deviation of solution of the BFKL equation in the low ^ region 
from the double logarithmic form (corresponding to (f) — 1). The curves 
correspond to different z values. At ,^ = 10, the curves increase with 
decreasing z from 2; = 1 to ^ = 0.1 in steps of 0.1. 

Fig. 4 The peculiarly normalised cross section for j*q — > Vq, calculated nu- 
merically by successively iterating the BFKL kernel. The dotted line 
corresponds to two gluon exchange, the dashed line incorporates the 
corrections predicted by the first iteration (of the BFKL kernel), the 
solid line is the prediction after the second iteration and the diamonds 
are the predictions after the third iteration. The dash-dot line is the 
asymptotic result obtained by computing eq.(9) numerically. No re- 
summation of the double logarithms has been performed, (a)^; = 0.05; 
(b) z = 0.35. 

Fig. 5 As in fig. 4 except that the double logarithms have been resummed 
in the manner described in the text. 

Fig. 6 As in fig.5, but for (a) z = 0.2; (b) z = 0.5; (c) z = 0.8. 
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